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Abstract. Isolated HAEBE stars are believed to represent an intermediate stage of objects between young stellar 
objects surrounded by massive, optically thick, gaseous and dusty disks and Vega like stars surrounded by debris 
disks. The star AB Auriga is already known for being surrounded by an intermediate-stage dust disk emitting 
a fairly large infrared and (sub-)millimetric excess. Until now, the outer disk structure has only been resolved at 
millimeter wavelengths and at optical wavelength coronographic imaging. We have obtained 20 /im images which 
show an unexpected ellipse-shaped disk structure in emission at a distance of about 260 AU from the central star. 
Large azimuthal asymmetries in brightness can be noticed and the center of the ellipse does not coincide with the 
star. A simple, pure geometrical model based on an emission ring of uniform surface brightness, but having an 
intrinsic eccentricity succeeds in fitting the observations. These observations give for the first time clues on a very 
peculiar structure of pre-main-sequence disk geometry, i.e. a non uniform increase in the disk thickness unlike 
the common usual sketch of a disk with a constant flaring angle. They provide also valuable informations on the 
disk inclination as well as its dust composition; at such a large distance from the star, only transient heating of 
very small particles can explain such a bright ring of emission at mid-infrared wavelengths. Finally, the increase 
of thickness inferred by the model could be caused by disk instabilities; the intrinsic eccentricity of the structure 
might be a clue to the presence of a massive body undetected yet. 
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1. Introduction 

Herbig Ae/Be (HAEBE) stars represent a class of interme- 
diate mass, pr e-main-sequen ce (PMS) stars, first described 
as a group bv lHerbid l)l960(l . The circumstellar (CS) disks 
found around these stars are believed to be the sites of 
on-going planet formation. By studying the characteris- 
tics and evolution of the CS disk and its dust composi- 
tion, valuable insights can be obtained into the processes 
leading to the formation of planets, and put constraints on 
disk and planet formation models. Infrared spectroscopy 
obtained with the Infrared Space Observatory (ISO) has 
given us insight into the dust composition of sample 
of isolated HAEBE systems (e.g. iBouwman et al.ll2000al 
[2001.: .Meeus et all Eo Ol). While these spectra reveal a 
rich mineralogy, no direct information concerning the spa- 
tial distribution of the different dust species can be in- 
ferred from the ISO data. Most studies so far have used 
the available spectral energy distributions (SEDs) to put 
constraints on the spatial distribution of the CS mate- 
rial. Models for passively heated disks surrounding PMS 
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stars are successful in reproducing the ISO s pectra (e.g. 
iDullemond et alJl200lt IPominik et al.ll20n3a^ ■ but these 
mode ls cannot be uniquely co nstrained from SED fitting 
alone iBouwman et al.l l(20fl0bl) . for this spatially resolved 
imaging, as presented in this paper, is required. 

Among isolated HAEBE systems, the disk around AB 
Aurigae is one of the most intere sting and studied. Its 
star has a probable age of 2 Myr l|van den Ancker et alJ 
119971^ . indicating that, according to current planet for- 
mation theories, planet building could still being ongo- 
ing in this system. ISO spectra show strong PAH emis- 
sion bands, and emission from silicates, and carbonaceous 
dust g rains ijvan den Ancker et al.ll200ol : IBouwman et alJ 
l2000al) . Though the inferred grain sizes are differing from 
interstellar grains, the dust around AB Aur seems to be 
relatively unprocessed, indicating an evolutionary young 
system. This seems also to be confirmed by combined 
ISO-SWS and sub-millimeter observations of H2 and CO 
rotational lines, d emonstrating th at the disk still has a 
large gas content (iThi et alJl200l|) . Its disk has been re- 
solved in the millimeter range by iMannings fc Sargend 
showing a structure consistent with a Keplerian 
disk. Its surrounding nebulae and outer disk structure 
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Fig. 1. Image of the AB Aurigae disk at 20.5 /im. The left figure shows the deconvolved image, with a pixel scale 
of 0.3"/pixel. Clearly visible is a resolved central emission region surrounded by a ring like structure. The panels on 
the right show the normalized intensity profiles along a cut through the CAMIRAS images of 2 reference stars and 
AB Aur. The upper right panel shows the comparison between two observations of reference stars, demonstrating the 
stability of the PSF. The panel on the lower right shows a comparison of AB Aur with the PSF, clearly showing that 
the central emission is extended and that a ring-like emission structure is also detected in non deconvolved data. 



were also studied in the vi sible range using broad-band 
coronographic observations <)Gradv et al showing a 

disk with s piral-shaped structures. Near-IR interferomet - 
ric (Millan-Gabet et all Il999l 1200 it lEisner et al.ll2003|) . 
have resolved the inner parts of the disk, showing it to 
be consistent with a passive disk with an inner hole, seen 
at a low (< 40°) inc lination angle. Recent mid-IR imaging 
I Chen fc ,Tural2003|) . has also resolved the inner structure, 
showing it to originate from thermal emission from dust 
grains heated by the stellar radiation field near the central 
star. Here we present mid-IR imaging, not only resolving 
the thermal emission from dust close to the central star, 
but also an emission structure at the outer parts of the 
disk. In the next sections we will describe these observa- 
tions and propose a pure geometrical model for explain- 
ing our results, unlike the paradigm according to which 
gravitational interactions with massive bodies are inferred 
when observing brightness asymmetries in dusty disks. 

2. Data collection and Reduction 

The observations w ere performed using the CEA mid-IR 
camera CAMIRAS l|Lagage et al.lllQQ'iji . equipped with a 
Boeing 128x128 pixels BIB detector sensitive up to a wave- 
length of « 28 fj,m. AB Aurigae was observed from the 
CFH 3.6m telescope on 1999 August 1 and between 2000 



March 16 and March 20. During these two runs, seeing and 
weather conditions - humidity and amount of atmospheric 
precipitable water - were particularly favorable and ex- 
tremely stable in time. We spent a total integration time 
on AB Aurigae of 2h split into several nights. Each dataset 
was reduced independently in order to avoid any erroneous 
conclusion due for instance to some corrupted dataset; it 
permits also to evaluate error bars on our results. The ori- 
entation of the array on the sky was carefully determined 
at the start of each observing run. The pixel size was 0.29". 
In the 20 /im window, we used a filter centered at 20.5 /im 
and with a bandpass (FWHM) A A = 1.11 fim. This fil- 
ter is free of any important atmospheric line contribution. 
The source was always observed at an airmass of less than 
1.3. Standard chopping and nodding techniques were ap- 
plied with a chopping throw of 16" North and nodding 
amplitude of 20" West. A shift-and-add procedure was ap- 
plied to each final cube of images using a correlation based 
method with a re-sampling factor of 4:1. In order to get 
the best spatial resolution in the dataset, the nodding di- 
rection was perpendicular to the chopping direction, re- 
sulting in an image containing four times the image of the 
source. The four sources were extracted and co-added us- 
ing a cross-correlation method. Finally, each independent 
dataset was deconvolved separately before co-addition to 
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Fig. 2. The structure of the emission ring as an function of azimuth, compared with our best-fit model. Plotted are, 
from left to right, top to bottom, with the solid lines, the ring radius, its width (fwhm), and its intensity profile as a 
function of azimuthal angle (from the west direction, going counter clock wise) . Over-plotted with dashed lines are the 
modeled ring radius, width and intensity profile. For the three solid line plots (data), the profiles have been smoothed 
by a 5 points (number of angle points being 36) boxcar in order to smooth high frequencies variations. The error bars 
indicate the formal 1 a error. Also shown is the best fit model image (bottom right) with overlaid the rotation axis 
(link to the disc inclination) and the major axis of the intrinsic ellipse described in sectional 



produce the final image. The peak signal-to-noise ratio in 
the various datasets ranges from 225 to 255. The standard 
star a Tau was frequently monitored for further data pho- 
tometric calibration and for PSF measurements. Because 
final CO- added images are strongly limited in spatial reso- 
lution by the seeing and above all by the 3.6m telescope 
diffraction-pattern (FWHM of 1.4" at 20.5 /im), the use 
of deconvolution techniques is mandatory in order to re- 
cover the best spatial resolution in MIR images. We used 
the Multiscale Maximu m Entropy Method developed by 
IPantin fc Starckl l)l996|) based upon wavelet analysis based 
on the concept of multiscale information. Shannon theo- 
rem prescription allows to pursue the deconvolution down 
to twice the pixel size, in this case down to « 0.6", and 
iterations are stopped according to the residual map (i.e. 
the data minus the re-convolved solution), which proper- 
ties should be consistent with noise characteristics. Each 
epoch dataset was processed (including image deconvolu- 



tion) independently in order to check the consistency of 
the results. 

3. Results and discussion 

3.1. Emission structures and photometry 

After deconvolution, the 20.5 /im image of AB Aurigae 
shows two main structures, as can be seen in the left panel 
of Fig.^ The most striking one is a narrow ring-like struc- 
ture slightly open (i.e. the flux is decreased by an order 
of magnitude) in the SW direction. It has a total flux of 
2.7±0.3 Jy (taking into account deconvolution errors) and 
is located at an average distance of 1.8 "(260 AU, assum- 
ing a distance of 144 pc) and with a typical width FWHM 
of 0.45 "(65 AU). The second structure, and brightest, 
with a total flux of 48.7±5 Jy, is concentrated in the close 
vicinity of the star. By comparing AB Auriga; and PSF 
profiles we find that the central emission is also resolved 
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Table 1. Best fit model parameters of AB Aur. Listed are 
the model parameters defining the system orientation, the 
ring emission and central emission region. 



Parameter Value ± la error 

system orientation! : 



e 


65 ° 


± 


2.0 ° 


i 


-9 ° 


± 


1.0 ° 


ring parameters: 






e 


0.13 


± 


0.01 


a 


1.95 " 


± 


0.01" 


$ 


129 ° 


± 


5.0 ° 


AR 


0.12 " 


± 


0.0034 " 


Ho 


0.57 " 


± 


0.038 " 


AH 


0.163 " 


± 


0.012 " 


Iring 


1.98 Jy arcsec"^ 


± 


0.017 Jy arcsec 


parameters 


central emission: 






Pccntcr 


0.42" 


± 


0.0044 " 


Iccntcr 


38 Jy arcsec"^ 


± 


1.4 Jy arcsec"^ 



jangles are given in degrees counter clockwise with re- 
spect to the west direction (right direction) 



with a FWHM (diameter) of 3 pixels or 0.9 "(130 AU), 
probably due to thermally emitting dust grains in the in- 
ner regions of the circumstellar disk (note that the star 
emission (w 4 mJy) is negligible at that wavelength). 

3.2. Interpretation: disk structure and the nature of 
the emitting grain population. 

The inner resolved structure, in our schematic view of the 
disk thermal emission, corresponds to thermal emission of 
hot grains at the inner parts of the disk, between 10 AU to 
about 65 AU from the star, where the temperature of the 
dust grains is such that they thermally emit efficiently at 
around 20 /iin . Wc find a similar structure as reported by 
IChen fc Jiiral |2003,) . which can be explained by thermal 
emission from either a disk or an envelope. We will dis- 
cuss and model this emission extensively in a forthcoming 
paper. 

Here we will concentrate on the intriguing ring-like 
structure at about 260 AU from the star. First of all, 
"standard" thermal emission from "big grains" is dis- 
carded at such a large distance from the star, be- 
cause the stellar radiation would heat them to a tem- 
perature around 100 K or less. These grains would 
then produce a flux in the range 30 to 60 fim in- 
consistent with the spectrum measured by ISO/SWS. 
Therefore, we must infer transient heating of very 
small (like PAH particles for instance) by UV radi- 
ation, that relax through emission of narrows bands. 
PAH bands were already found in the SWS spectrum of 
AB Au rigae ijvan den Ancker et alJl2000t iBouwman et all 
2000a). Inde e d, acc ording to ISchutte et alJ (|l993^ and 
Draine fc Lil l)200lh . large PAH molecules (or small 



grains), consisting of a few thousand carbon atoms, can 
efficiently emit infrared radiation with a "plateau" -like 
spectral shape around 20 /im producing then the observed 
emission of w 3 Jy at 260 AU from the star; derived rela- 
tive band strength are also consistent with PAH emission 
models cited above. 

The fact that the geometrical emission is arranged in a 
ring puts strong constraints on the disk geometry. In the 
following we call apparent or projected ellipse the ellipse 
directly seen on the 20 iim image. The most constrain- 
ing feature of this ring-like elliptical structure is that it 
is off-centered along the major axis of the projected el- 
lipse. In principle, inclining a circular structure leads to a 
projected elliptical structure that can be oflF-centered but 
always along the apparent minor axis. Also the systematic 
brightness variation along the ring is not symmetric with 
respect to the minor axis of the apparent ellipse meaning 
that either this brightness variation is intrinsic to the ring 
or that the orientation of the inclination axis is not the 
same as the major axis of the ellipse. This leads unavoid- 
able to the conclusion that we are looking at an intrin- 
sically elliptical structure seen under a small inclination 
angle. 

To determine the main characteristics of the ring struc- 
ture as seen in the CAMIRAS image of AB Aurigae, we 
constructed a simple, geometrical model based on the de- 
convolved image. Since the inferred very small grains need 
stellar UV or visible photons to be excited, their emission 
most likely originates from the disk surface layer, i.e. the 
disk photosphere, being directly illuminated by the central 
star. To mimic this emission from the disk surface layer, we 
used a uniform brightness ring, with an eccentricity e, semi 
major axis a, and orientation $, radial width AR, scale- 
height Hq, vertical width AH, and surface brightness Iring- 
The ring model is positioned such that the central emis- 
sion is in the focus of the ellipse. The central emission is 
modeled with a uniform brightness disk with radius Reenter 
and surface brightness Iccntcr- Further, the disk model has 
a overall inclination i and orientation O. This "infinite 
resolution" mode is then convolved with a cr=0.25 gaus- 
sian filter to take into account the limited resolution of 
the deconvolved image (0.6"). Our best-fit model (using 
the 2D fitt ing procedure MP FIT2DFUN of IDL package 
written bv lMarkwar'd^l)2002|) ) is shown in Fig. Inland the 
resulting fit parameters are listed in tabled We checked 
also that the overall inverse modeling is correct by verify- 
ing that our solution corresponds also to an minimum of 

in raw data (non deconvolved) space. A sketch of the 
derived disk structure is shown in Fig.|3| The appearance 
of an emission ring can be linked to the disk geometry as 
follows: In a flaring disk model the amount of radiation 
intercepted and thus emitted is proportional to the angle 
at which the stellar light impinges onto the disk surface. 
This angle is a factor of 15 (with the parameters found 
in the case of AB Aurigae) larger at the ring surface com- 
pared to a uniform flaring (wedge) disk. This implies that 
a factor of 15 more light is intercepted and consequently 
emitted at the disk surface, lifting it above the detection 
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Fig. 3. Schematic representation of the disk structure in 
the AB Aur system. Drawn is the derived geometry of 
the disk surface. Indicated in the figure are the model 
parameters defining this geometry as hsted in Table ^ 
Rring is the distance to the star and varies between a(l — e) 
and a(f + e). 



limit of the CAMIRAS instrument. This large flaring an- 
gle of the disk surface can also explain the gap in the SW 
direction, as the disk surface becomes self shadowed seen 
under a small inclination of —20°. This kind of "puffed- 
up" structure could be for instance the re sult of shadow- 
ing inst abilities as described in paper by iDominik et alJ 
Other possibilities to explain a sudden increase 
of brightness at a distance of 260 AU, would be: 



— dust sizes segregation l|Takeuchi fc Linll2002j) but one 
must infer a very peculiar structure of the gaseous disk 
to reproduce the observations. 

— PAH photo-dissociation, but it is unlik ely to occur 
to the inferred large PAH grains (e.g. lAllain et al.l 

— an hollow spherical shell and its limb brightening, but 
the off-centering can be hardly explained. 

Remains to explain the origin of the intrinsic eccentric- 
ity of the ring-like structure assumed in the model de- 
scribed above. Some massive perturber on an eccentric 
orbit might produce the detected gravitational deforma- 
tions of the ring. A giant planet is relatively unlikely to 
be found one at such a high distance; a not yet detected 
brown dwarf could be one possible explanation. 

4. Conclusions 

We have presented in this paper a thermal infrared im- 
age at 20.5 of the AB Aurigae dusty disk. The decon- 
volved image shows an inner resolved structure contain- 
ing the majority (95 %) of the thermal flux at this wave- 
length and an unexpected ring-like ellipse-shaped struc- 
ture. This ring is located at an average distance of 260 
AU from the star, but its most interesting features lie in 
large scale asymmetries and a center that is offset with 
respect to the star. When modeling the ring structure 



using a purely geometrical model (in which asymmetries 
are essentially produced by inclination and shadowing ef- 
fects in a flared disk), we end with the conclusion that we 
are in presence of an intrinsically ellipse-shaped ring 
that could trace the presence of puffed-up matter from a 
flared dust disk containing transiently heated particles. 
Self-shadowing instabilities perturbating the disk vertical 
thickness could produce such a ring. However, the origin 
of its intrinsic eccentricity remains unclear; massive bod- 
ies could be acting in gravitationally structuring the disk. 
Further detailed modelings still need to be investigated to 
fully understand the physics of this intriguing structure; 
further observations, including higher resolution imaging 
in several PAH bands and spatially resolved mid-IR spec- 
troscopy, are needed to better understand this object. 
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